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q ' Abstract. Space-based gravitational- wave interferometers such as LISA will be sensi- 

tive to the inspiral of stellar mass compact objects into black holes with masses in the 
£^ ' range of roughly 10 5 solar masses to (a few) 10 7 solar masses. During the last year of 

inspiral, the compact body spends several hundred thousand orbits spiraling from sev- 
eral Schwarzschild radii to the last stable orbit. The gravitational waves emitted from 
these orbits probe the strong-field region of the black hole spacetime and can make 
possible high precision tests and measurements of the black hole's properties. Measur- 
ing such waves will require a good theoretical understanding of the waves' properties, 
which in turn requires a good understanding of strong-field radiation reaction and of 
. properties of the black hole's astrophysical environment which could complicate wave- 

form generation. In these proceedings, I review estimates of the rate at which such 
inspirals occur in the universe, and discuss what is being done and what must be done 
further in order to calculate the inspiral waveform. 
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One of the most exciting sources that should be measured by the space-based 
gravitational-wave interferometer LISA (the Laser Interferometer Space Antenna 
[1]) is the inspiral of a "small" (1 — 10 M©) compact body into a massive (10 5-7 M & ) 
black hole. Such massive black holes reside at the cores of galaxies; the smaller 
compact bodies will become bound to the hole and spiral into them after undergoing 
interactions with stars and other objects in the environment of the galactic center. 
The measurement of gravitational waves from such inspirals will make possible very 
high precision tests of general relativity, and probe the nature of the galactic core's 
environment. 

To set the stage for understanding why these inspirals are such interesting objects, 
consider the following estimates: the orbital energy of a small body in an equatorial, 
prograde orbit of a Kerr black hole is 



£° rb = /x }~K +qv3 , , (i) 



where v = yM/r and q = a/M. (I use units in which G = c = 1 throughout.) The 
orbital frequency of the small body is 



ft = 



M 1 ' 2 



(2) 



r 3 / 2 + aM 1 / 2 • 



Radiation reaction carries orbital energy away from the system, causing the orbit 
to shrink. Eventually it shrinks enough that the body reaches the innermost stable 
circular orbit (ISCO). Orbits inside this radius are dynamically unstable; further 
radiative evolution tends to push the body into the black hole. 

Post-Newtonian theory allows us estimate the rate at which the system loses 
energy as a power series in the quantity u = (Mft) 1 / 3 (which is roughly the orbital 
speed) and the black hole's spin a. Reference [2] gives the energy loss in such a 
post-Newtonian expansion: 



The prefactor — 32/5(/x/M) 2 -u 10 is the result one gets applying the quadrupole for- 
mula to a binary described with Newtonian gravity; the function /schw.(w) is a 
(rather high-order) correction appropriate for zero-spin black holes, and f sp i n (a,u) 
is a correction incorporating information about the hole's spin. (Note that this 
formula is only appropriate for fj, <C M: it does not incorporate any finite mass 
ratio corrections.) 

Equations (l)-(3) can be used to estimate the time it takes for a small body to 
spiral from radius r\ to r 2 , and the number of gravitational- wave cycles it emits in 
that time: 



(On the last line, I have assumed that the bulk of the radiation comes out in the 
quadrupole m = 2 mode, so that / gw = 2ft/27r.) Consider now the inspiral of a 5 M 
body into a rapidly spinning (a ~ M) 10 6 M black hole. The small body spirals 
from a radius of 8M (in Boyer-Lindquist coordinates) to the ISCO in one year, 
emitting around 5 x 10 5 gravitational-wave cycles as it does so. The gravitational- 
waves that it emits lie in the frequency band 3 x Hz < / < 3 x 10~ 2 Hz — the 
band to which LISA is most sensitive. Indeed, careful analyses of the detectability 
of the signal by LISA [3] indicate that such an inspiral should be detected out to a 
distance of roughly 1 Gigaparsec with amplitude signal-to-noise ratio of around 10 
to 100 (depending on factors such as the mass of the small body, the mass of the 
black hole, and the black hole's spin). The fact that such a large number of cycles 
are emitted indicates that details of the waveform can in principle be measured 
to very high precision. Detection of extreme mass-ratio inspirals by LISA offers 
the possibility of very high precision measurements of the characteristics of extreme 
strong-field regions of spacetime. 
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The high precision tests that LISA will be able to make should allow us to directly 
map the characteristics of the massive body's spacetime metric and confirm that 
they in fact exhibit the Kerr metric. Most likely the way that this will be done 
will be to measure the multipole moments of the massive body. Fintan Ryan [4] 
has shown that the gravitational waves which are emitted as a small body spirals 
into a massive compact object contain a "map" of the massive body's spacetime. 
By measuring the gravitational waves and decoding the map, one learns the mass 
and current multipole moments which characterize the massive body. All multipole 
moments of Kerr black holes are parameterized by the holes' mass M and spin a: 

M t +iSi = M (ia) 1 . (6) 

For Kerr black holes, knowledge of the moments M = M and Si = aM deter- 
mines all higher moments. This is one way of stating the "no-hair" theorem: The 
macroscopic properties of a black hole are entirely determined by its mass and spin. 
(I neglect the astrophysically uninteresting possibility of charged holes.) By mea- 
suring gravitational waves from extreme mass ratio inspiral and thereby mapping 
the massive body's spacetime, LISA will test the no-hair theorem for black holes, 
determining whether the massive body has multipole moments characteristic of the 
Kerr metric, or whether the body is something more exotic, such as a boson star. 

The waves emitted by extreme mass ratio inspiral are thus likely to be directly 
measureable by LISA, and are likely to be extremely interesting. One might next 
wonder whether they occur often enough to be interesting. This question has 
been examined in some detail by Martin Rees and Steinn Sigurdsson [5,6]. They 
consider the scatter of stellar mass black holes in the central density cusp of galaxies 
into tightly bound orbits of the galaxies' central black hole. Occasionally, such a 
scattering event will put the stellar mass hole into an orbit which is so tightly 
bound that its future dynamics are driven by gravitational-wave emission, and it 
becomes an interesting source for LISA. They find that the rate of such events is 
likely to lie in the range 

1 event ^ „ ^ 1 event 
<TZ< 

year Gpc 3 ~ ~ month Gpc 3 

Obviously, there are large uncertainties in this calculation. The low mass end of the 
massive black hole population (which is most relevant to LISA observations) is not 
as well constrained as the population of very massive black holes (M > 1O 8 M ), 
and there are uncertainties in the rate at which stellar mass black holes are "fed" 
into the central hole to produce extreme mass ratio binaries. However, the lower 
end of the rate estimate (7) is based on very conservative estimates. We may rather 
robustly estimate that the rate measured by LISA will be several events per year 
out to a Gigaparsec [7]. 

The waves that LISA will measure from these sources will come from orbits 
that are rather eccentric and inclined with respect to the black hole's equatorial 
plane [8]. To best interpret the measured waves (and, indeed, in order to improve 
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the odds of seeing the waves at all in the detector's noisy data stream), it will be 
necessary to have some theoretical modeling of the orbit and the waves that it emits 
as gravitational radiation reduces the orbit's energy and angular momentum. We 
expect that the radius, inclination angle, and eccentricity of the orbit will change as 
radiative backreaction drives the system's evolution. Some means of understanding 
these changes, in detail, is needed in order to model the waves' evolution accurately. 

Before discussing work in radiation reaction, it is necessary to take a sanity check. 
Relativity theorists often work in a very idealized universe: their extreme mass ratio 
binary is likely to consist of a big black hole, a small body, and gravitational waves. 
In the real astrophysical world, there will be complications to this pretty (but 
highly idealized) picture. One should worry whether the complications render the 
relativity theorist's modeling invalid. 

Perhaps the most important such complication arises from interaction between 
the small inspiraling body and material accreting onto the massive black hole. 
Recently, Narayan has analyzed this interaction and concluded that, in almost all 
cases, accretion induced drag is unlikely to significantly influence extreme mass 
ratio inspiral [9]. This conclusion is based on the fact that in the majority of 
cases, the rate at which the central black hole accretes gas from its environment 
is rather low (several orders of magnitude less than the Eddington rate [10]). For 
these "normal" galaxies, much evidence [9,10] suggests that the gas accretes via 
an advection dominated accretion flow (ADAF). Narayan's calculation [9] shows 
that the timescale for ADAF drag to change the orbit's characteristics (e.g., the 
orbital angular momentum) is many (9 to 16) orders of magnitude longer than 
the timescale for radiation reaction to change the orbit's characteristics. Thus, 
the relativity theorist's idealized view of an extreme mass ratio binary is probably 
quite accurate: radiation reaction is likely the most important factor driving the 
evolution of extreme mass ratio binaries. 

When the mass ratio is extreme, one can analyze the spacetime of the binary using 
a perturbative expansion: the spacetime metric can be written as a "background" 
from the central object (which I will assume from now on is a Kerr black hole), 
plus a perturbation due to the inspiraling body: 

g a/3 = g^(M,a) + h al3 (tx). (8) 

The evolution of the perturbation h a p(n) should then describe the dynamical evo- 
lution of the system. To linear order in the mass ratio fi/M (which should be 
adequate for extreme mass ratios), this evolution can be described using perturba- 
tion techniques, such as the Teukolsky equation 1 [11]. 

When the mass ratio is extreme, the effects of radiation reaction are gentle enough 
that the system's evolution is adiabatic: the radiation reaction timescale is much 
longer than the orbital timescale. At any given moment, the trajectory of the 
small body is very nearly a geodesic, parameterized by the three constants of Kerr 



^ The Teukolsky equation actually describes the evolution of a curvature quantity related to the 
perturbation. 



orbital motion: the energy E, the (^-component of) angular momentum L z , and the 
"Carter constant" Q. Gravitational-wave emission causes these three constants to 
change on the radiation reaction timescale. In this adiabatic limit, the evolution of 
the system can be understood in terms of the evolution of the quantities (E, L z , Q). 

It is well known that gravitational waves carry energy and angular momentum. 
One might think that they carry "Carter constant" as well, and that therefore one 
might be able to deduce the effects of radiation reaction by measuring the flux of 
radiation at infinity and down the event horizon. By measuring the amount of E, 
L z , and Q carried in the flux one should be able to deduce how much E, L z , and 
Q are lost from the orbit. This would then allow one to figure out orbits of Kerr 
black holes radiatively evolve. 

This approach does not work. One can deduce the change in the orbit's E and 
L z by examining radiation fluxes, but one cannot so deduce the change in Q: 

8E or bi t = — 5 -Eradiated , 
^^z, orbit radiated j 

^Qorbit 7^ ~~ Eradiated • (9) 

The change 5Q turns out to depend explicitly on the local radiation reaction force, 
= dp>*/dT, which the small body experiences due to radiative backreaction (see 
[14]). The properties of this force (and programs to calculate it) are described 
elsewhere in this volume [12]. Here, it is sufficient to note that an understanding 
of the radiation reaction force for gravitational radiation reaction lies some time in 
the future, so that we cannot yet evolve generic Kerr black hole orbits. 

There are special cases where evolution of the Carter constant is not such a nasty 
impediment. One case is the evolution of equatorial orbits. Equatorial orbits have 
Q — 0; and, one can easily show that an orbit which starts off equatorial remains 
equatorial. In this case, one need only evolve the energy and angular momentum. 
The local radiation reaction force is not needed in this case. Another case is the 
evolution of circular, non-equatorial orbits. (For non-zero spin, "circular" means 
"constant Boyer-Lindquist coordinate radius".) Such orbits have recently been 
shown to remain circular under adiabatic radiation reaction [13]. Thus, in an 
adiabatic evolution, a system which is initially circular and inclined will remain 
circular and inclined: the system evolves through a sequence of orbits changing 
only its radius and inclination angle. By imposing "circular goes to circular", one 
can write down a simple rule relating the change of the Carter constant to the flux 
of energy and angular momentum: Q = Q(E,L Z ). 

Recently, I have examined the evolution of these circular orbits under adiabatic 
radiation reaction, using a flux-measuring formalism based on the Teukolsky equa- 
tion. The formalism and results are presented at length in [14]. Some highlights of 
the results are presented in Figure 1. 

Consider the left panel of Figure 1. The horizontal axis is orbital radius r; the 
vertical axis is inclination angle l. This figure shows the direction, in (r, t) phase 
space, in which radiation reaction tends to push the orbit. This particular analysis 




FIGURE 1. The left panel illustrates the effect of radiation reaction on circular, inclined orbits. 
It is for a black hole with a = 0.8M. The tail of each arrow represents a particular orbit with 
radius r and inclination angle t. The arrow indicates the direction in which radiation reaction 
drives the orbit; the size of the arrow indicates how quickly it is so driven. The right panel is a 
snapshot of the gravitational waveform emitted during an inspiral into a hole with a = 0.95M. 
At this moment, the small body is spiraling through r = 7M, i = 62.43°. 



is for a black hole with a = 0.8M. Notice that the orbits are in the strong-field of 
the hole; the dotted line indicates the maximum inclination angle which the orbit 
can have and remain stable. Orbits tilted beyond this line are dynamically unstable 
to small perturbations and plunge into the hole. The tail of each arrow represents a 
particular orbit. The direction of the arrow gives the direction in which that orbit 
tends to evolve due to gravitational-wave emission; the arrow's length indicates 
the relative rate of evolution. In all cases, the direction of the arrow is such that 
the inclination angle increases: radiation reaction tends to make tilted orbits more 
inclined. This is exactly what one would have predicted by extrapolating from 
post-Newtonian theory [15]. The rate at which this inclination angle increases is 
rather slow — the aspect of the arrows in Figure 1 is nearly flat. Indeed, the 
value of di/dt in this strong-field regime is roughly 3 times smaller than what post- 
Newtonian theory predicts 2 . An interesting feature of Figure 1 is the very long 
arrow at r = 7M, i ~ 120°. This orbit lies extremely close to the marginally stable 
orbit: it is at i = 119.194°; the marginally stable orbit is at i = 119.670°. This 
orbit is barely dynamically stable, so a small push has drastic effects. 

The right panel of Figure 1 shows a portion of the gravitational waveform emitted 
during inspiral. The central black hole in this case has a = 0.95M; it lies at 



2 ) Recent analyses are showing that in the extreme strong-field of rapidly rotating holes, the incli- 
nation angle changes rather more dramatically, and in the opposite direction: radiation reaction 
tends to decrease the inclination angle. This work is in progress [16]. 



luminosity distance D from the detector. The waveform here is shown as the small 
body passes through r = 7M, i = 62.43°. Note the low-frequency modulation 
of both polarizations. This is due to the frame dragging induced precession of 
the orbital plane — Lense-Thirring precession. Note also the many sharp, short- 
timescale features present in the two polarizations. When the spin is high, many 
harmonics of the small body's fundamental orbital frequencies contribute to the 
gravitational waveform. This leads to a rather complicated structure; the energy 
spectrum corresponding to these waveforms extends to rather high frequencies. 
Accurate measurement of such complicated waveforms will be quite a challenge. 
However, the payoff is likely to be immense. 
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